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Abstract 



Measurements of structure functions and of the hadronic final state in 
deep inelastic scattering at HERA are presented. The results comprise the 
extraction of parton densities, measurements of the strong coupling, and the 
search for novel QCD effects in the new kinematic regime at HERA. 
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1 Introduction 



1.1 Overview 

Data from the electrons-proton collider HERA at DESY in Hamburg provide new 
insights into the structure of the proton. Electrons with E e = 27 GeV collide with 
protons of E p = 820 GeV , giving rise to a total ep centre of mass (CM) energy of 
^fs = ^AE e E p pa 300 GeV . The high energy allows to probe the proton in hitherto 
unexplored kinematic regions. Two complementary sources of information can be 
identified: the measurement of the inclusive ep — *> e'H deep inelastic scattering (DIS) 
cross section, where H stands for any hadronic system, and measurements of the 
hadronic system H itself. The cross section measurements are expressed in structure 
functions, which in turn can be interpreted in terms of parton (i.e. quark and gluon) 
densities in the proton. In simple words, how often the probe hits something inside 
the proton determines its content. 

The hadronic final state H on the other hand can be viewed as the materialization 
of the quantum fluctuations inside the proton, allowing in principle a more direct, 
but due to confinement also a more complicated, access to the dynamics inside the 
proton. In addition, due to the large phase space, the hadronic final state serves as 
a QCD laboratory with tunable initial kinematic conditions. 

1.2 Kinematics 

The kinematics of the basic scattering process in Fig. [I] can be characterized by 
any set of two Lorentz-invariants out of Q 2 ,x,y,W, which are built from the 4- 
momentum transfer q = k — k! mediated by the virtual boson, and by the 4- 
momentum p of the incoming proton. These are: 

Q 2 = —q 2 , which gives the resolving power of the probe with wavelength A = 
1/Q (we set Ti = c = 1); 

x = Q 2 /(2pq), the Bjorken scaling variable, which can be interpreted as the 
momentum fraction of the proton which is carried by the struck quark; 

y = Q 2 I (xs) , the transferred energy fraction from the electron to the proton in 
the proton rest frame; and 

W 2 = Q 2 (l — x)/x + m 2 k, sy — Q 2 , the invariant mass squared of the hadronic 
system H . 

The large CM energy give access to kinematic regions both at very small x 
and at large Q 2 (Fig. g). The HERA data cover roughly Q 2 = 0.2 - 10 4 GeV 2 , 
x = 10~ 5 - 10- 1 and W = 40 - 300 GeV . 



2 here the generic name electron includes the positron 
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Figure 1: a) Basic diagram for DIS in 0(a®) (quark parton model - QPM). b) As an example, 
a DIS event with the scattered electron and two well separated jets detected in the HI detector is 
shown. The proton remnant leaves mostly undetected in the +z direction. 

1.3 Detectors 

Two large multi-purpose colliding beam detectors, HI and ZEUS, are installed 
around the collision regions (Fig. [l|b). They consist mainly of inner tracking de- 
vices, followed by calorimetry and outer muon detection systems. They allow to 
measure both the scattered electron and the hadronic final state, apart from the 
proton remnant, which largely escapes detection in the beam pipe. The scattered 
electron serves as a tag of DIS events, distinguishing them from an otherwise over- 
whelming background from photoproduction (quasi real photons with Q 2 « 0) and 
beam-gas reactions. The kinematics can be determined either from the electron 
alone, or from the measured hadronic system alone, or from a combination of both, 
permitting important systematic cross checks. 



2 Structure functions 

The fundamental measurement in DIS concerns the cross section for ep — > e'H as 
a function of the kinematic variables. The quark parton model (QPM) offers a 
physical picture: the scattering takes place via a virtual photon which is radiated 
off the scattering electron, and which couples to a pointlike constituent inside the 
proton, that is a quark or antiquark. The cross section is then proportional to the 
quark density inside the proton. The general expression for the differential cross 
section ep — > e'H is 



dV 



4na z 



dxdQ 2 xQ 4 



l- y +?L\.F 2 (x,Q 2 ) 



j-Fl(x,Q 2 ) 



(1) 



where Z exchange has been neglected (it is a 1% correction for Q 2 « 1000 GeV 2 ). 
The structure function F 2 can be interpreted (in the "DIS" scheme) in terms of the 
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Figure 2: Coverage in the kinematic plane (x, Q 2 ) of various DIS experiments. 

quark and antiquark densities, and q~i, and their couplings to the photon, i.e. their 
charges e qi : 

F 2 (x,Q 2 ) = xJ2e 2 q Mx,Q 2 ) + q l (x,Q% (2) 

i 

where the sum runs over all quark flavours. The longitudinal structure function Fi 
vanishes in zero'th order a s , and will be discussed later. 

In the simple quark parton model the proton consists just of 3 valence quarks. 
Their distribution functions in fractional proton momentum x , xq(x), would peak 
at ~ 1/3 and tend towards zero for x — > 0, 1. In a static model of the proton, they 
would not depend on Q 2 . It then follows that F 2 should not depend on Q 2 , just 
on x (Bjorken scaling). 

When QCD is "turned on" the quarks may radiate (and absorb) gluons, which 
in turn may split into quark - antiquark pairs or gluon pairs. More and more 
of these fluctuations can be resolved with increasingly shorter wavelength of the 
photonic probe, A = l/Q. With Q 2 increasing, we have a depletion of quarks at 
large x , and a corresponding accumulation at lower x. In addition, "sea quarks" 
from g qq splittings populate small x. In fact, at small x it is the gluon content 
with distribution function g(x, Q 2 ) which governs the the proton and gives rise to 
the DIS cross section via qq pairs. 

The dynamic features of DIS have been calculated in various approximations 
of QCD. One obtains equations for the Q 2 dependence of parton densities (and 
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structure functions) like 



dq(x,Q 2 ) _ a s (Q 2 ) fidz 
ding 2 ~ 2tt 



q(z,Q 2 )P qq (^)+g(z,Q 2 )P w ( 



(3) 



(and similarly for the gluon densities), the famous DGLAP (Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi) equations [0. They involve the calculable Altarelli-Parisi 
splitting functions P{j(z), which give the probability of parton branchings q — > qg, 
g — > gg and g —>■ qq, where the daughter parton i carries a fraction 1 — z of the 
mother's (j) momentum. The coupled integro-differential equations for the quark 
and gluon densities can be solved, allowing to calculate them for any value of Q 2 , 
once they are known at a particular value Qq. 

For large Q 2 and for small x where the proton is dominated by gluons one obtains 
in the "double leading log" (DLL) approximation (lnln(Q 2 / Qq) ^> 1, ln(l/x) 3> 1) 
the formula 



xg{x,Q 2 ) w xg(x,Q 2 )exp 



144 ln^Mln(l/,) 



\ 25 ln(Q 2 /A 2 



(4) 



At small x a fast rise of the gluon density is predicted. That is, xg increases faster 
than (ln-) A , but slower than (^) A for any power A, the growth depending on the 
"evolution length" from Ql to Q 2 . 



2.1 The structure function F 2 

ZEUS and HI have measured F 2 in a completely new kinematic domain compared 
with fixed target experiments, most notably to much smaller values of x. The 
HERA data || || (Fig. ^) exhibit a steep rise of F 2 towards small x , which flattens 
at smaller Q 2 . They signal increasing parton densities at smaller x. 

In Fig. [| F 2 is shown as a function of Q 2 for fixed x values. At x ~ 0.1 scaling 
is observed - F 2 does not depend on Q 2 . At x > 0.1 F 2 decreases with Q 2 due to 
parton splittings, the products of which are found then at smaller x . Therefore 
F 2 increases with Q 2 for x < 0.1 in accord with the above qualitative discussion. 

From a DGLAP evolution of pre-HERA data this sharp rise could not be pre- 
dicted a priori, because input distributions at small x were not available. It was 
known though that asymptotically for Q 2 — > oo the small x behaviour is given 
by the DLL formula eq. |4j. The interesting question is whether the HERA data 
are consistent with DGLAP evolution, or whether new effects have to be taken 
into account, which one may expect at very small x. The DGLAP evolution re- 
sums all leading terms ~ (a s \n(Q 2 / Ql)) n , but neglects terms ~ (a s ln(l/x)) n 
in the perturbation series. The latter may become important at very small x 
(a s ln(l/x) ^> 1) but moderate Q 2 , and are treated in the BFKL (Balitsky-Fadin- 
Kuraev-Lipatov) equation. The BFKL equation predicts an even more dramatic 
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growth (the "Lipatov growth") of the gluon density than the DLL expression at 
small x, F 2 ~ xg(x) ~ x~ A , with A ~ 0.5 in leading order (LO). 

It turns out, however, that the data can be fit perfectly well with parton densi- 
ties which obey the next-to-leading-order (NLO) DGLAP evolution equations (see 
Fig. ^|) []§, U . Standard QCD evolution works over many orders of magnitude in both 
x and Q 2 \ HI obtains values between 0.2 and 0.4 for the exponent in F 2 ~ (l/x) A , 
increasing with Q 2 . In fact, both the x and the Q 2 dependence of F 2 can be at- 
tributed predominantly to the DLL formula, which can be displayed nicely with a 
suitable variable transformation ("double asymptotic scaling" @). 

It is remarkable that the F 2 data can be described by evolving flat or valence- 
like input quark and gluon distributions from a very low scale, Q 2 , = 0.35 GeV 2 
up in Q 2 H . The steep rise with decreasing x is achieved by the long evolution 
length in Q 2 from the very low scale Ql = 0.35 GeV 2 (see eq. |J). The success 
of the GRV prediction came as a surprise for many, as perturbative QCD should 
only be applicable for Q 2 not too close to Aq CD « 0.2 GeV , because otherwise 
a s (Q 2 ) = 127t/ ((33— 2nj) ln(<5 2 /AQ CD )) diverges. It is quite satisfying that the data 
can be described with parton densities following standard QCD evolution. However, 
the goal remains to calculate the measured magnitude of the growth (l/x) A from 
QCD, rather than tuning it with the starting point Qq of the DLL evolution. 



2.2 QCD analysis of F 2 — the gluon density and a s 

From the DGLAP eq. [| it is clear that the scaling violations of F 2 depend on both 
a s and the gluon density. In fact, for x < 0.01 and in lowest order one can derive 
the approximate formula [ l2| 



dF 2 (x/2,Q 2 ) 10 a s (Q 



dlnQ 2 27 7T 



xg(x,Q ), (5) 



because at small x the proton is dominated by gluons, and the scaling violations arise 
from quark pair creation from gluons. The full NLO QCD analysis now employed 
at HERA is of course more involved. Fig. ||]a shows the gluon density x ■ g(x,Q 2 ) 
extracted from the NLO QCD fit to the F 2 data. It increases sharply towards small 
x, the rise being more pronounced at larger Q 2 . 

Clearly, the density cannot increase forever; eventually saturation effects will 



have to be taken into account [13, |14|. When gluons of transverse size ~ 1/Q fill 
up the whole transverse area offered by the proton, they will start to overlap and 
recombine. This would be a novel and very interesting situation indeed: high parton 
density, but Q 2 large enough for a s (Q 2 ) to be small! Given the size of the proton 



lfm, the critical condition for saturation effects to turn on can be estimated HI4" 



as x cr i t (7(x cr it, Q 2 ) ~ yj^z ~ 25 ■ This value is by far not reached by the 

measured gluon density (Fig. ||a). It could be however that saturation does not set 
in uniformly over the proton's transverse area, but starts locally in so-called hot 
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Figure 3: The structure functions F%(x, Q 2 ) as a function of x for different Q 2 bins. Shown are 
data from HI g, ZEUS f§, BCDMS §, NMC @ and E665 (llj. A NLO DGLAP QCD fit to 
the data with Q 2 > 5 GeV 2 is shown as full line. 



spots [p~5|| . In this would be larger. The data however do not require any 

saturation correction. 

From an analysis of the scaling violations di^ /dlnQ 2 ~ a s , or equivalently, of 
a QCD fit to the F2 data, the strong coupling constant can be determined. From an 
analysis of the 1993 data a s {mz ) = 0.120±0.005(exp.)±0.009(theor.) was obtained 
[M . This is being updated for the higher statistics data now available p 



2.3 The longitudinal structure function Fl 

The structure function F 2 can also be expressed in terms of the cross sections ot and 
a 1, for the absorption of transversely and longitudinally polarized virtual photons^] 
on protons: a^t = (Jl + c"r , namely 

„ _ q\i-x) q 2 



47r 2 a Q 2 + 4m 2 x 2 



a 1 p Sri 
"tot ~ 



4vr 2 







(a L + (T T ), 



(6) 



3 we use the Hand convention 111] for the definition of the virtual photon flux 
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Figure 4: The structure functions F2(x, Q 2 ) as a function of Q 2 for different x bins. For better 
visibility, a function C(x) — 0.6(i — 0.4), with i the x bin number (i = 1 for x = 0.32), is added to 
F2 in the plot. Shown are the HERA and fixed target data, as well as a NLO QCD fit {§. 



where the small x approximation has been applied. Similarly, 

r q\i-x) _ g 2 



4ir 2 a 



4ir 2 a 



■ °l ■ 



(7) 



Longitudinal photons have helicity and can exist only virtually. In the QPM, 
helicity conservation at the electromagnetic vertex yields the Callan-Gross relation, 
Fl = 0, for scattering on quarks with spin 1/2. This does not hold when the quarks 
acquire transverse momenta from QCD radiation. Instead, QCD yields |L9] 



4:71 Jx Z 6 



+ 8j2el(l-- z )-g(z,Q 2 
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exposing the dependence of Fl on the strong coupling and the gluon density, which 
dominates over the first term at small x. In fact Fl(x,Q 2 ) = 0.7 • [4a s /(37r)] • 
xg(2.5x, Q 2 ) is not a bad approximation for x < 10~ 3 pQ| . 

The extraction of F 2 from the cross section measurement (eq. [1|) so far had to 
make an assumption for Fl, because there existed no data in the HERA regime. At 
large y, y pa 0.7, this is a 10% correction. The argument can be turned around, and 
Fl can be extracted at large y from a measurement of the cross section, assuming 
that F 2 follows a QCD evolution and can be extrapolated from measurements at 
smaller y. This procedure has been performed by HI |2T| (Fig. |5|b). The result 
Fl = 0.54 ± 0.03 ± 0.22 excludes the extreme possibilities Fl =F 2 and Fl =0 and 
implies R = (Tl/o't = Fl/(F 2 — Fl) ~ 0.5, since F 2 pa 1.5. This is self consistent 
with the gluon density extracted from the HI QCD fit. 
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Figure 5: a) The gluon density xg(x, Q 2 ) as extracted from a QCD analysis of the HI F2 data, b) 
The longitudinal structure functions Fl ■ The error bars represent statistical and systematic point- 
to-point errors. The error band gives an additional global systematic error. The QCD expectation 
based on the GRV M parton densities is also shown (solid line, unlabelled). 



2.4 The very low Q region 

The ZEUS collaboration had for the 1994 run a new beam pipe calorimeter installed 
to cover small electron scattering angles. They can now access Q 2 values below 
1 GeV 2 which had previously been the realm of fixed target experiments, see Fig. 



[ 22| . Will perturbative QCD still hold at such a low scale? The perturbatively 
evolved parton densities of GRV can describe the data down to Q 2 ~ 1.5 GeV 2 . 
Below Q 2 ^0.4 GeV 2 the GRV curves turn over and become valence like, by far 
failing to account for the data. A non-perturbative model by Donnachie-Landshoff 
|23| , which is based on Regge phenomenology and reproduces the total photoproduc- 
tion cross section (Q 2 pa 0) gives a good description of the data in this regime, but 
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fails above Q 2 « 2 GeV 2 . The transition from perturbative to non-perturbative 
QCD appears to be rather fast. 



• ZEUS 95 BPC(preliminary) A ZEUS 94 □ E665 DL GRV94 
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Figure 6: The structure function F 2 measured at small Q 2 by ZEUS p2[ . Overlayed are the 
expectations from the Donnachie-Landshoff (^J Regge model (DL), and from the perturbatively 
evolved parton densities from Gliick, Reya, Vogt || (GRV). 



3 The hadronic final state 

New insights into both perturbative QCD and confinement can be gained by study- 
ing the hadronic final state, thereby providing complementary information to the 
totally inclusive structure function measurements. For example, footprints of BFKL 
dynamics are being searched for in especially devised observables (Et flows, forward 
jets, high p T particles). F 2 may be too inclusive to be sensitive to these effects, be- 
cause all possible paths in the parton evolution are integrated over. Differences 
between DGLAP and BFKL evolution can be expected for the hadronic final state, 
because in DGLAP the phase space for parton radiation, basically given by W, is 
restricted by strong ordering of the transverse momenta kx of subsequently emitted 
partons, which is not the case for BFKL (Fig .0a). 

There exists also a possibility to discover new QCD effects, for example spec- 
tacular QCD instanton events. At a more "down to earth" level, the large phase 
space for hard QCD radiation (W up to 300 GeV) can be exploited for the extrac- 
tion of the strong coupling constant from jet rates and other final state obervables. 
For both cases, the searches for new effects and precision measurements, a good 
understanding of the final state dynamics in parton showers and hadronization is a 
prerequisite for meaningful results 

Apart from the laboratory frame, the hadronic centre of mass system (CMS) 
and the Breit frame are used in the analyses. The Breit frame is defined by the 
condition that the virtual photon does not transfer energy, only momentum. In 
the QPM picture the scattering quark would thus just reverse the direction of its 
momentum of magnitude Q/2. The CMS is defined as the centre of mass system 
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Figure 7: a) A ladder graph for parton evolution. The selection of forward jets in DIS events 
is illustrated, b) The cross section for forward jets as a function of x. Overlayed are QCD model 
predictions for hadron jets, and BFKL and Born graph calculations for partons. The jets are 
selected with x ]ct > 0.035, 0.5 < k T j e t 2 /Q 2 < 2 and fc Tjct > 3.5 GeV . 



of the incoming proton and the virtual boson, i.e. the CMS of the hadronic final 
state with invariant mass W. In both systems the hemisphere defined by the virtual 
photon direction is referred to as the current region, the other (containing the proton 
remnant) as the target region. The CMS current and target systems are back to back 
with momentum W/2 each. Longitudinal and transverse quantities are calculated 
with respect to the boson direction. 



3.1 QCD models 

The complexity of the hadronic final state makes analytical QCD calculations of- 
ten very difficult - there exist but a few. Indispensable tools for the study of the 
hadronic final state are therefore Monte Carlo models which simulate the under- 
lying dynamics. They incorporate QCD evolution and parton radiation in differ- 
ent approximations and utilize phenomenological models for the non-perturbative 
hadronization phase. The MEPS model (Matrix Element plus Parton Shower), an 
option of the LEPTO generator P^j , incorporates the QCD matrix elements up to 



first order, with additional soft emissions generated by adding leading log parton 
showers. In the colour dipole model (CDM, generator ARIADNE) (ESI ^| radia- 
tion stems from a chain of independently radiating dipoles formed by the colour 
charges. The HERWIG model |27| is also based on leading log parton showers, with 



additional matrix element corrections [28] . The CDM description of gluon emission 



is similar to that of BFKL evolution to the extent that the gluons emitted by the 
dipoles do not obey strong ordering in hr |[29|| . In the MEPS and HERWIG models 



the partons emitted in the cascade are strongly ordered in kx, because they are 
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based on leading log DGLAP parton showers. 



3.2 Transverse energy flow 



It can be shown |K| that DGLAP evolution corresponds to the evaluation of ladder 
graphs (Fig. [?[) in which subsequently emitted gluons are strongly ordered in trans- 
verse momenta hp , Ql <C £>r \ ■■■for j ^ ■■■Q 2 ■ m the BFKL evolution on the 
other hand gluon emission is not restricted by strong ordering, the emitted gluons 
rather follow a kind of random walk in kp space. As a consequence, one expects 
more gluon radiation and hence more transverse energy Et emitted in the rapidity 



region between the remnant and the current systems [31]. The HERA data |3 



are compatible with the predicted magnitude |H[ and the characteristic increase to- 



wards small x from the BFKL calculation (Fig. ^). However, measured are hadrons, 
but the calculations do not include hadronization! It turns out that DGLAP models 
with large hadronization corrections to the partonic Et are also able to describe the 



data [34, pq|, precluding strong conclusions on BFKL dynamics from the Et data. 




Figure 8: a) The average hadronic Et vs. x, measured in —0.5 < rj < 0.5 in the CMS from 
HI [H and ZEUS [f| (both preliminary). The BFKL calculation at the parton level is also 
shown, b) The charged particle pt spectra measured in 0.5 < i] < 1.5 in the CMS by HI for high 
and for low x at Q 2 14 GeV , and the predictions from the colour dipole model (ARIADNE), 
LEPTO and HERWIG. The plots are normalized to the total number of events N. 



11 



3.3 Inclusive transverse momentum spectra 

Another method to probe the low x parton dynamics is to look for high p? particles 



36fl . They signal the presence of hard gluons, while hadronization can only produce 
limited pt • Perturbative parton radiation can thus be disentangled from hadroniza- 
tion effects, which pose a problem for the interpretation of the Et data. HI finds 
at low x and at central rapidity significantly harder pt spectra than expected from 
Monte Carlo calculations based upon LO matrix element and leading log DGLAP 
parton showers (Fig. ^o). The colour dipole model, in which parton radiation is 
more abundant, gives a good description of the data. QCD calculations and BFKL 
Monte Carlo predictions for these spectra are eagerly awaited to find out whether 
the enhanced parton radiation can indeed be attributed to BFKL evolution, or to 
something else, perhaps even more interesting. 

3.4 Forward jets 



A classic signature for BFKL dynamics |37|, |38| is the production of "forward jets" 
with Xj e t = Ejet/Ep, the ratio of jet energy and proton beam energy, as large as 
possible, and with transverse momentum hp^ close to Q in order to reduce the 
phase space for the ftp ordered DGLAP evolution (s. Fig. [7]). An enhanced rate of 
events with such jets is thus expected in the BFKL scheme J37|. The experimental 
difficulty is to detect these "forward" jets which are close to the beam hole in the 
proton direction. The rate of forward jets measured by HI |39| (Fig. |^b) increases 
with falling x. This is expected from BFKL calculations, in contrast to calculations 



without the BFKL ladder [f40| . The behaviour of the data is better represented by 



the CDM than by the DGLAP representative, the MEPS model. However, the cross 



sections are calculated 40 for partons, while experiments measure hadron jets. This 



gap has to be bridged from both sides to allow a strictly valid comparison. 
3.5 Inst ant ons 

The standard model contains processes which cannot be described by perturbation 
theory, and which violate classical conservation laws like baryon and lepton number 
in the case of the electroweak sector and chirality for the strong interaction jLL]. Such 



anomalous processes are induced by instantons f42fl . At HERA, QCD instantons may 



lead to observable effects in the hadronic final state in DIS [H3L The instanton 
should decay isotropically into a high multiplicity state, consisting of gluons and all 
quark flavours (in each event!) which are kinematically allowed. Due to the isotropic 
decay, one expects a densely populated region in rapidity, other than the current 
jet, which is isotropic in azimuth. The presence of strangeness and charm provides 
an additional signature (Fig. ||a). HI does not see an excess of K° production over 
the prediction from standard QCD models (Fig. ||) |45|]. At most a few percent 



admixture of instanton events is allowed to normal DIS, corresponding to an upper 
limit for the instanton cross section for x > 0.001 of 0.95 nb. In the future, more 
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elaborate search strategies and larger luminosities offer the chance for a fundamental 
discovery at HERA f46fl ! 




Figure 9: a) Instanton induced process, and the K° yield as a function of laboratory pseudora- 
pidity 7^ The proton direction is to the right. Shown in b) are the predictions from standard QCD 
models, and from an instanton model (QCDINS). In c) the HI data 0] are overlayed with the 
prediction of a standard QCD model plus a 0.6% instanton contribution, and with the maximally 
allowed~jjraction of instanton gvents (QCDINS). 



3.6 The running coupling constant ot s from jet rates 



The processes contributing to jet production in DIS up to first order in a s are shown 
in Figs. [1,10. The QPM process results in a so-called "1+1" jet topology, while the 
QCDC and BGF processes give "2+1" jet events, where the "+1" refers to the 
unobserved nemnant jet. Since the rate of 2+1 jet production is oc a s , a s (Q) can 
be measured |47], |48| , and the scale dependence can be studied. The ZEUS data 
of 1994 are shown in (Fig. [Top. The Q dependence is as expected for a running 



a s . A fit yields a s (m z ) = 0.117 ± 0.005(stat.)l^(syst.) ± 0.007(theor.), to be 
compared with the world average, a s {mz) = 0.117 ± 0.005. It can be expected that 
the theoretical uncertainties can be reduced by studying different jet algorithms j|9 
and by obtaining a better understanding of the hadronic final state dynamics. 



3.7 Thrust in the Breit frame 

In the Breit frame in- and outgoing quark have equal but opposite sign momenta 
Q/2 (QPM picture), and in e + e~annihilation the outgoing quark and antiquark have 
equal but opposite momenta \fs/2 = Q/2. Due to this similarity it is interesting to 
compare measurements in the Breit current hemisphere in DIS with e + e~data. The 
current hemisphere of the DIS Breit frame is equivalent to one hemisphere of the 
e + e~ — > qq reaction. DIS experiments have the advantage that they cover a large 
span in Q 2 in a single experiment. 

Infrared safe event shape variables have in the past been a very useful way to 
measure ot s , because it was possible to calculate them in QCD. One example is the 
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Figure 10: Determination of the strong coupling constant. Shown are the boson-gluon fusion 
(BGF) and the QCD Compton (QCDC) graphs, giving rise to 2+1 jet events, and the measured 
a s as a function of Q from the ZEUS jet rates p8|. 



thrust T, where T is the normalized sum of all longitudinal particle momenta with 
respect to the thrust axis ny: 

rp E \Pi ■ n T \ 

1 = max ■ 



40 50 60 70 f 



Ht is varied to maximize the thrust. T = 1 when all particles are collinear, and 
T = 0.5 for an isotropic distribution. The measured thrust value in the Breit 
current hemisphere [|50| increases with increasing Q (i.e. increasing energy of the 



scattered quark) - the current jet becomes more collimated (Fig. [11]). The data agree 
well with what has been measured in e + e~experiments. The advantage at HERA 
is evident: the evolution with Q can be studied in a single experiment. In fact the 
data can be fit with a QCD ansatz (1 — T) — c\ • a s (Q) + c 2 • a s 2 {Q) + c 3 • 1/Q, 
where the first terms are perturbative up to NLO, and the last term parametrizes 
all higher orders and non-perturbative hadronization effects with a power correction 



51]. The HERA data will allow to check the hypothesized universality of the power 



correction, and to extract a s (Q) once the NLO coefficient has been calculated. 



3.8 Scaling violations of charged particle spectra 

The transition from partons to hadrons can be described with fragmentation func- 
tions D(z, Q), which give the probability to find a hadron carrying the momentum 
fraction z of the original parton's momentum. These functions are universal, and 
exhibit scaling in lowest order, i.e they are independent of Q. In DIS one can mea- 
sure the momentum spectra in the scaled variable x p = p/(Q/2)), where p is the 
hadron momentum and Q/2 the maximal possible momentum in the Breit frame. 
When QCD radiation is switched on, the original quark may radiate a gluon, and 
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e 




Figure 11: QCD analysis in the Breit frame (see sketch), a) The average 1 - T (T=thrust) 
measured in the Breit frame as a function of Q. The data are compared to e + e _ data and to the 
LEPTO program, both for hadrons and for partons. A QCD fit to the data is also shown, b) The 
normalized charged particle cross section in bins of x p = 2p/Q as a function of Q 2 [p^| . 



instead of one hard parton to fragment there are now two softer ones, which result 
in a softer x p spectrum. These scaling violations increase with Q, as the phase space 
for QCD radiation increases. The preliminary ZEUS data pi on x p (Fig. |11|) indeed 
show that with increasing Q there are less hadrons with large x p , and more hadrons 
with small x p . These scaling violations are analogous to the scaling violations of 
Fo , and can be used in a similar fashion to extract a, 151 



4 Conclusions 

HERA has already delivered a wealth of information about the structure of the 
proton, in particular about partons which carry a small fraction of the proton's mo- 
mentum. Standard QCD evolution gives a good description of the inclusive structure 
function measurements. Novel QCD effects are being searched for in the hadronic 
final state with some interesting signals emerging, the interpretation of which is un- 
der discussion. In addition, the strong coupling constant is being determined using 
a variety of hadronic final state observables. These analysis profit from the fact 
that the scale at which the coupling constant is being determined can be varied in 
a single experiment at HERA. 
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